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ABSTRACT. 

Shock  Models  Arising  from  Processes  with  Stationary,  Independent,  Non¬ 
negative  Increments 
Harry  Joe  and  Frank  Proschan 

Let  H(t)  be  the  life  distribution  of  a  device  subject  to  shocks 
governed  by  an  integer-valued  stochastic  process  with  stationary,  in¬ 
dependent,  nonnegative  increments.  H(t)  is  a  function  of  the  probabi¬ 
lities  ^of  surviving  the  first  k  shocks,  k  •  1,  2,  ...  .  We  show  that 
H(t)  inherits  various  aging  properties  (2FR,  I FRA,  NBU)  of  the  discrete 

Jl 

survival  function  F.  .  Analogous  results  hold  in  the  continuous  case 
where  H(t)  is  the  life  distribution  of  a  device  subject  to  wear  accord¬ 
ing  to  a  wear  process  with  stationary,  independent  nonnegative  incre¬ 
ments.  In  the  cumulative  damage  model  for  the  wear^rocess  or  for  the 
shock  process,  H  is  I  FRA  if  the  process  enjoys  a  Tf^property. 


✓:i  tAi  ice: 
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1.  Introduction  and  Summary. 

In  Esary,  Marshall,  and  Proschan  (1973)  (hereafter  referred  to  as 
EMP  (1973)),  the  life  distribution  H(t)  of  a  device  subject  to  shocks 
governed  by  a  Poisson  process  is  considered  as  a  function  of  the  proba¬ 
bility  of  surviving  the  first  k  shocks,  k  =  1,  2,  ...  .  EMP(1973) 

show  that  the  continuous  life  distribution  H(t)  inherits  various  aging 
properties  of  the  discrete  survival  function  such  as  IFR,  I FRA,  NBU, 
etc. 

In  Section  2,  we  extend  some  of  these  results  to  the  case  in  which 
shocks  occur  according  to  an  integer-valued  stochastic  process  (N(t), 
t  a  0}  with  stationary,  independent,  nonnegative  increments  (i.e.,  a 
generalized  Poisson  process) .  In  Section  3,  similar  results  are  ob¬ 
tained  in  the  continuous  case  of  a  wear  process  (W(t) ,  t  a  0)  with  wear 
occurring  according  to  a  stochastic  process  with  stationary, independent, 
nonnegative  increments.  In  Section  4,  we  present  some  examples  of  well 
known  generalized  Poisson  processes.  Finally,  in  Section  5,  we  study 
cumulative  damage  models  for  the  shock  process  (N(t),  t  a  0)  and  the 
wear  process  (W(t),  t  a  0).  In  the  cumulative  damage  model  for  the 
shock  process,  the  amounts  of  damage  caused  by  the  shocks  are  indepen¬ 
dently  distributed  according  to  a  common  distribution,  and  the  device 
fails  when  the  total  accumulated  damage  exceeds  a  specified  threshold. 

In  this  model,  the  resulting  probabilities  P^  enjoy  the  discrete  I FRA 
property,  and  thus  the  life  distribution  H(t)  of  the  device  is  IFRA. 

In  the  cumulative  damage  model  for  the  wear  process,  the  amount  of  da¬ 
mage  caused  by  a  fixed  amount  of  wear  always  has  the  same  distribution, 
and,  as  before,  the  device  fails  when  the  total  accumulated  damage 

! 

i 


exceeds  •  specified  threshold.  In  this  nodel,  the  resulting  probe- 
bilities  Fu  of  surviving  the  first  u  units  of  wear  enjoy  the  I FRA 
property  (under  one  assumption)  and  thus,  the  life  distribution  H(t) 
is  again  1FRA. 
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2.  Generali led  Poisson  Process 

Let  a  device  be  subject  to  shocks  that  are  governed  by  a  genera¬ 
lized  Poisson  process.  A  generalized  Poisson  process  is  an  integer¬ 
valued  stochastic  process  (N(t),  t  2  0}  with  stationary,  independent, 
nonnegative  increaents  (see  Parzen,  1962,  p.  126) .  Suppose  the  pro¬ 
bability  of  the  device  surviving  k  shocks  is  7^,  where  1  ■  Fq  a  7j  2 
P i  *  *"•  It  follows  that  the  survival  function  FT(t)  of  the  device  is 
given  by: 

ff(t)  »  Jo  Pjj  PCN(t)  -  k).  (2.1) 

2.1.  Definition.  A  process  (N(t),  t  2  0}  is  said  to  be  TP2  if 
PCN(t)  ■  kj  is  TP2  in  t  e  [0,  •)  and  k  e  (0,  1,  2,  ...}.  (See  Karlin, 
1968,  for  the  definition  of  TP2.) 

The  following  theorem  extends  a  basic  result  of  EMP(1973) . 

2.2.  Theorem,  (i)  Let  7.  be  discrete  NBU  (F.  s  77  for  m,  k  ■ 

— —  k  K+m  K  in 

0,  1,  2,  ...).  Then  H  given  in  (2.1)  is  NBU.  (ii)  Let  7^  be  discrete 

IFR  (7.  is  log  concave  in  k  «  0,  1,  2,  ...)  and  let  (N(t),  t  2  0}  be 

1/k 

TP2.  Then  H  is  IFR.  (iii)  Let  7^  be  discrete  IFRA,  i.e.,  7^  is  de¬ 
creasing  in  k  «  0,  1,  2,  ...  .  Let  (N(t),  t  2  0}  be  TP2.  Then  H  is 
IFRA. 

Proof.  Let  t,  s  2  0.  Then 

k 

!T(t  ♦  s)  «  £  7.  p[N(t  +  s)  -  k3  ■  Z  7.  I  P[N(t)  -  mJPCN(s)  -  k 

k«0  K  k-0  *  a»0 

-  PCN(t)  «  m]  A'*  PCN(s)  -  k  -  mJ 
•  ■  Jo  p£N(t)  -  I  7.  PCN(s)  -  kj. 

m-0  kmQ  k+m 
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The  second  equality  follows  from  the  fact  that  the  increments  are 
stationary  and  independent. 

(i)  Since  F^  is  discrete  NBU,  then 

00  00 

H(t  ♦  s)  S  I  PCN(t)  -  m)  Pm  r  F.  P[N(s)  -  k] 
m»0  "  k«0  K 

•  H(t)  H(s), 

so  that  H  is  NBU. 

(ii)  Since  F^  is  discrete  IFR ,  then  F^  is  log  concave  in  k,  or  equiva¬ 
lently,  P^+B  is  TP2  in  k  and  -m.  Since  P[N(s)  *  k]  is  TP2  in  k  and  s,  it 

follows  from  the  Basic  Composition  Formula  (Karlin,  1968,  p.  17)  that 
00 

£  F.  P[N(s)  *  k]  is  TP2  in  s  and  -m.  Applying  the  Basic  Composition 
k*0  K  m 

Formula  once  more,  and  the  fact  that  P  [N(t)  =  m3  is  TP2  in  -m  and  -t,  we 
conclude  that  H(t+s)  is  TP2  in  s  and  -t.  Thus  H  is  log  concave  and  therefore 
IFR. 

(iii)  Let  QN^tj  (♦)  denote  the  probability  generating  function  of 
N(t).  Then  for  0  s  C  s  1, 

00 

1  (V  PWW  *  ^3  *  H(t)  -  (O 

k»0  k  N<tJ 

-  H(t)  -  CQN(1)(c)3t  -  H(t)  -  e_tC'l0g  Vl)003* 

The  second  equality  follows  from  the  relationship  between  the  probability 
generating  function  and  the  moment  generating  function.  Note  also  that 
f or  0  s  ?  s  1,P[N(1)  ■  03  <  *  1  »so  that 

Os-  log  QNrn(0  S  -log  P  CN(1)  -  03. 

N(1)  1/k 

Since  F^  is  discrete  IFRA,  then  F^  is  decreasing  in  k.  It  follows 

—  k 

that  P^  -  c  changes  sign  at  most  once,  and  if  one  sign  change  does  occur. 


it  occurs  from  *  to  -,  for  0  £  c  *  1.  Since  PC (N(t)  »  k]  is  TP2,  by 

•0t 

the  Variation  Diminishing  Property  (Karlin,  1968,  Chapter  S),  H(t)  -  e 
changes  sign  at  most  once,  and  if  once,  from  +  to  -,  for  0  <  0  £  -log  P[N(1)»0]. 
Note  also  that 

H(t)  2  P[N(t)  -  0]  -  QN(t)(0)  -  0j(1)(0)  «  e't{'l0g  PCN(1)  "  0J} 

It  follows  that  H(t)  2  e'0t  for  6  2  -log  P(N(1)  »  03.  Thus  H(t)  -  e"6t 
changes  sign  at  most  once,  and  if  once,  from  ♦  to  -,  for  6  >  o.  By  a  charac¬ 
terization  of  IFRA  distributions  (Barlow  and  Proschan,  1975,  p.  89),  H 
must  be  IFRA,  | | 

2.3  Remark.  In  a  dual  fashion,  similar  results  hold  for  the  dual 
classes  NWU,  DFR,  DFRA.  The  NBUE  and  DMRL  classes  and  their  duals  have 
not  been  successfully  treated. 


i 

I 

i 

} 

? 
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3.  Continuous  Wear  Processes 


In  this  section  we  obtain  analogous  results  for  continuous  wear 
processes.  Let  (W(t),  t  a  0}  be  a  stochastic  process  with  stationary, 
independent  nonnegative  increments  with  density  function  f^^C*)* 
Suppose  now  a  device  is  subject  to  wear  in  accordance  with  this  sto¬ 
chastic  process.  Let  Fu  be  the  probability  of  surviving  a  cumulative 
wear  of  u  units,  with  Pu  satisfying  1  =  i  Pu  i  Pu*  for  0  <  u  £  u". 
We  express  the  survival  function  IT(t)  of  the  device  as 


H(t)  *  Iq  Pu  fw(t)(u)  du*  (31) 

3.1  Definition.  The  process  (W(t),  t  i  0}  is  said  to  be  TPj  if 
fw(t) (°)  is  TP 2  in  t,  u  e  [0,  »). 

We  may  now  state  and  prove  the  continuous  analogue  of  Theorem  2.2. 

3.2  Theorem,  (i)  Let  Fu  be  NBU.  Then  H  given  in  (3.1)  is  NBU. 
(ii)  Let  Pu  be  IFR.  Let  (W(t),  t  *  0}  be  TPr  Then  H  is  IFR.  (iii) 
Let  ?u  be  IFRA  and  (W(t),  t  z  0}  be  Tp2*  Then  H  is  IFRA. 

Proof.  Let  t,  s  2:  0.  Then  H  (t  ♦  s)  *  Q  Pufw(t+s)(u)  du 


*  r0  puf»(o  *  f»(s) <u)  du 

•  to  Tufo  fw(s),u  -  x)  fw(t)(x)  dx  du 

■  to  •  x)  Fu  du  dx 

'  Q  f»(t)(X)  Q  Vs)00  du  dx- 

(i)  and  (ii)  now  follow  as  in  Theorem  2.1. 

(iii)  Let  Mw(tj(*)  be  the  moment  generating  function  of  W(t). 
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Then  for  0  <  c  <  j£(Pu  -  e*cu)fw(t) (u)du  -  H(t)  -  /”e*5ufw(t) (u)du 


-  H(t)  -  Mw(t)(-0  «  H(t)  -  tMw(1)(-C)]t  »  H(t)-  e 


-tC-iog  Mw(i) C-0 3 


For  0  <  z,  <  »,  we  have  0  <  £  1,  so  that  0  £  -log  <  •• 

Completing  the  proof  as  in  Theorem  2.2  (iii),  we  conclude  that  H  is  I FRA.  || 
3.3  Remark.  As  in  the  discrete  case,  dual  results  hold  in  the  con¬ 
tinuous  case  for  NWU,  DFR,  DFRA.  Results  have  not  yet  been  obtained  for 
the  NBUE,  DMRL  classes  and  their  duals. 
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4.  Examples. 

A  generalized.  Poisson  process  (Parzen,  pp.  126-127)  is  an  integer¬ 
valued  stochastic  process  (N(t),  t  s  0}  with  stationary,  independent, 
nonnegative  increments.  In  a  generalized  Poisson  process,  events 
can  occur  simultaneously.  A  generalized  Poisson  process  necessarily 

has  a  characteristic  function  of  the  form  4^^ (u)  =  e^^CuJ-l'J  for 

00  iku 

some  X  >  0  and  some  characteristic  function  <J»(u)  ^^o^k  e  .  This 

form  of  characteristic  function  for  N(t)  arises  when  times  of  possible  events 

follow  a  Poisson  process  with  intensity  X  and  the  number  of  events 

occurring  at  such  times  are  independently  distributed  with  probability 

atoms  {p^}. 

If  we  choose: 

.  iu  n 

(i)  $(u)  =  e^  '  \  the  Poisson  characteristic  function, or 

(ii)  $(u)  *  (p(l-qelu)  *}r,  the  negative  binomial  characteristic  function, 
or 

(iii)  4>(u)  *  (pelu  +  q)n,  the  binomial  characteristic  function, 

we  obtain  a  generalized  Poisson  7?2  process  with  corresponding  frequency 
function: 

(i)  PCN(t)  =  k]  f  Mil  or 

k!  j*0  j! 

(ii)  PCN(t)  «  k)  «  qke'Xt  \  (rj*k'1)  — ^P^ ,  or 

j“0  K  j! 

(iii)  P[N(t)  -  k]  -  pke‘Xt  \  ("j)  qnj"k. 

j-0  k  j ! 

An  example  of  a  stochastic  process  (W(t),  t  z  0}  with  stationary,  inde¬ 
pendent  nonnegative  increments  such  that  the  density  f^^j(u)  is  TP^  in  t, 
u  c  (0,  •»)  is  the  gamma  process: 

fw(t)(u)  *  Xotuot  *e  Xu/r(ot),  o  >  0. 


iO 

5.  Cumulative  Damage 


We  next  consider  a  model  in  which  wear  causes  damage,  damages 
cumulate,  and  the  device  fails  when  a  specified  cumulative  damage  thres¬ 
hold  is  exceeded.  IVe  assume  that  the  amount  of  damage  caused  by  a  fixed 
amount  of  wear  always  has  the  same  distribution.  Let  the  amount  of 
damage  caused  by  one  unit  of  wear  have  distribution  F,  infinitely  divi¬ 
sible  and  satisfying  F(0“)  *  0.  Let  the  damage  threshold  be  denoted  by  x. 


and  the  probability  of  surviving  u  units  of  wear  be  denoted  by  F  .  Then 
Pu  ■  F^(x),  defined  more  precisely  below.  We  shall  show  that  under  an 
appropriate  assumption,  PU^U  *  lF^U"*(x)]^u  is  decreasing  in  u  2  0  for 
fixed  x,  so  that  Pu  has  the  continuous  IFRA  property. 

Let  X  (•)  be  the  unique  comp lex- valued  function  of  t  with  X(0)  *  0 
that  is  continuous  at  0  satisfying  $(t)  *  e*^  (see  Chung,  p.  241),  where 
$  is  the  characteristic  function  corresponding  to  F.  Since  $  is  an  in¬ 
finitely  divisible  characteristic  function,  4»r ( » )  =  er^ is  a  characteristic 

Crl 

function  for  all  r  2  0  (see  Parzen,  p.  124).  Let  F  •'be  the  distribution 

I* 

function  with  characteristic  function  <(>  .  Then 


F(r)  t  p(s)  _  p(r+s)  r>  s  >  0 

C  r) 

5.1  Theorem.  Let  x  be  a  continuity  point  of  F  for  all  r  2  0. 
frl  l/r 

Then  [F  (x))  is  decreasing  in  r  2  0. 


Proof.  For  an  arbitrary  distribution  function  G  with  G(0“)  *  0, 

[G^ (x) is  decreasing  in  positive  integer  n  for  each  fixed  x  2  0 

(see  Barlow  and  Proschan,  1975,  p.  94),  where  G^  is  the  n-fold  convolution 

ml  ®2 

of  G.  Let  rj  <  r2  be  positive  rationals,  Tj  —  ,  r2  — ,  nij  <  n^,  m2, 

n  are  positive  integers.  Then  applying  the  referenced  theorem  with  G  ■  F^n^ , 
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we  have: 

{ (F ^ l^n) )  ("2 ^  (x)  } ly,®2  *  {(F^1/,n^)^ml^(x)>l^ml  or 

{F(m2/n)(x)}n/m2  s  {F(ml/n)(x)}n/ml  or 
{F(V(x)}1/r2  i{P^«)1/rl 

Therefore,  (F^Cx)}*^1  is  decreasing  in  r  over  the  positive  rationals 
for  each  fixed  x  £  0. 

f  rl 

Now  suppose  x  is  a  continuity  point  of  Fv  for  each  r  2  0.  Let  r 
and  s  be  positive  reals,  r  <  s.  Let  {rn>  be  a  sequence  of  rationals  such 

that  r  +  r  and  r  <  s  for  all  n;  let  {s  }  be  a  sequence  of  rationals  such 

m  n  n  <0  1/r.  <«*>  l/< 

that  s^  +  s.  Then  rn  £  s^  for  all  n  implies  that  Cf  "(xJJ1'1*1  z  [F  (x)J  '  " 

for  all  n.  By  the  Continuity  Theorem, 

F^fx)  -»>  F(r)(x)  and  F(Sn)(x)  -*•  F(s)(x).  Thus  CF(rn) (x) ]1/rn 

-*•  CF(r)(x)]1/r  and  [F(^(x)  ]1/sn  -►  [F(s) (x) J1/s,  with  CF(r)(x)J1/r  > 

[F^(x)]^S.  Thus  we  conclude  {F^(x)}*^r  is  decreasing  in  r  2  0.  || 

Using  the  IFRA  property  of  F  and  of  F.,  we  obtain  the  following  theorems. 

5.2  Theorem:  Let  H(t)  =  /*  F^  (x)fw^t^  (u)  du  represent  the  survival  pro¬ 

bability  in  the  cumulative  damage  model  for  the  wear  process  (w(t),  t  2  0}. 
Suppose  x  is  a  continuity  point  of  F^  for  all  u  2  0  and  suppose  (w(t),  t  2  0} 
is  TPj.  Then  H  is  IFRA.  ..  . 

Proof :  This  follows  from  theorem  5.1  and  Theorem  3.2  iii. 

00 

5.3  Theorem.  Let  TT(t)  *  £  G^(x)P(N(t)  =  k)  represent  the  survival 

k=0 

probability  in  the  cumulative  damage  model  for  the  shock  process  (N(t),  t  2  0}, 
where  G  is  the  distribution  function  for  damage  due  to  a  single  shock  (G(0“)  ■  0) 
and  G^  is  the  k-fold  convolution  of  G.  Suppose  (N(t),  t  2  0}  is  Tp2* 

Then  H  is  IFRA. 

Proof.  This  follows  from  Theorem  2.2(iii)  and  the  fact  that  CG^(x))1^  is 
decreasing  in  k  for  each  x  2  0. 
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